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SECTIOM 1 
INTRODUCTION 


Laem 


In predicting rocket nozzle performance and analyzing Arc Plasma 

Generator (APG) data, the heat and mass transfer coefficient must be defined. 

In the past, the heat transfer coefficient has been calculated within 
Aerotherm by some boundary layer integral codes such as ARGEIBL (Reference 1) : 
or BLIMP (Reference 2). Both of these computer codes assume that the boundary 

layer develops on a smouth wall. It has been shown (see Reference 3) that 

surface roughness developed in rocket nozzle and APG environments can cause f 
suostantial heat transfer augmentation. To account for this and other 
boundary layer effects, the Momentum/Energy Integral Technique (MEIT) computer i 
program has been developed. . i 


ee, 


The MEIT methodology was originally developed to predict nosetip i 
ablation on reentry vehicles. This methodology has been validated for : 
reentry environment and {fs currently an important part of the ASC code i 
(Reference 4). The methodology employed by MEIT solves both the integral 
momentum and energy equations. In addition, it utilizes heat transfer and 
skin friction laws based on the most recent roughness heating and friction 
data. 
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MEIT was developed from ASC by extracting al] the subroutines that are 
related to the boundary layer calculation. Any nonapplicable subroutines, 
unnecessary variables and common block variables were deleted and a driver 
routine was written to read in all the required input. To tailor it for 
rocket nozzle applications, the input and output formats were modified so 
that they are compatible with the usage requirements of rocket nozzle designers. 


MEIT was verified by comparing soluttons to solutions generated by 
ARGEIBL and BLIMP. Solutions from the TBL code (Reference 5) which is 
developed by Pratt and Whitney Aircraft Company, were also used. Both 
"smooth and rough wall solutions were considered in the MEIT checkout. Like x 
ARGEIBL and other integral techniques, the heat transfer coefficient 
generated by MEIT must be decreased by 25 percent. 


In Section 2, a brief description of the MEIT methodology is presented. 
This is followed by the verification of MEIT in Section 3. Sectfon 4 provides 
the input and output formats and Section 5 presents two sample problems. A 
listing of the program is given in Section 6. 
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SECTION 2 
MEIT METHODOLOGY 


The MEIT methodology was criginally developed to predict nosetip 
ablation in reentry vehicles. A detailed description of the methodology is 
given in Reference 6. Only the essential information will be presented in 
Section 2. . 


The MEIT methodology solves both the boundary layer integral momentum 
and energy equation.. The required input are surface shape, boundary layer 
edge conditions , boundary layer gas properties and wall conditions (see 
Section 4). To solve these two equations, the local shape factor , recovery 
factor, Stanton number and friction coefficient must be defined. The effects 
of surface roughness, transpiration, acceleration, an¢ boundary layer 
properties are taken into account in terms of influence coefficients. These 
influence coefficients are included in the formulation of both the local 
Stanton number and friction coefficient. The solution procedure is carried 
out by an implicit finite difference scheme. Although MEIT fs designed 
primarily for rocket nozzle environment, which consists of turbulent flow; 
both laminar and transitional flow situations are also included. 


In Section 2.1, the basic equations are presented. This is followed 


by the formulation of local shape factor, recovery factor, Stanton number, 
and friction coefficient in Section 2.2 The solution procedure is described 
fn Section 2.3. 


2.1 BASIC EQUATIONS 


The two boundary layer integral equations solved by MEIT are: 


Integral momentum equation 


Cc (pv)_u 
gd ¢. af, we, HO dp 
ro .Ue ds (ri Qu,0) aes Oeue Oue 


Integral energy equation 


“ . iy 
moN), vite 
vy) Petelty eg = FY) 


J d 
PU ho ST as (PPctelhe g - WYO) = C it 
rel, hee hy, ds ee’ t,e w h hee -h 


bnere the momentum and energy thicknesses are respectively: 


The boundary layer shape factor, H, is defined ac: 


ne o& 


where 5*, the displacement thickness is given by: 


(1) 


hi) 


(4) 


5) 


(2) 
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while the recovery enthalpy is given by 


Dr 


i 


h, = he t Fa (8) 


where F is the recovery factor. 
The heat transfer rate and skin friction are related to the Stanton 


number and friction coefficient respectively by: 


(9) 


eR (10) 


In order to facilitate the solutions of equations (1) and (2), besides 
the required input of surface shape, bouncary laver edge conditions, boundary 
layer gas properties, and wall conditions, the local shape factor, recovery 
factor, Stanton number and friction coefficient have to be formulated. 


2.2 SHAPE FACTOR, RECOVERY FACTOR, STANTON NUMBER AND FRICTION COEFFICIENT 


The shape and recovery factors are evaluated in MEIT by the follow- 
ing relations. 


For laminar flow: 


I 
WwW 
i <2 3.029 + - 0.614 (11) 
: L (12) 
Fo Pr? 
For turbulent flow: 
t 
H, = 2.285 (1 + 3.2e7") r - 0.96 (13) 
W743 = 
Be 2 iPr (14) 


where n, the turbulent velocity profile exponent, is given by: 


0.37 + 2n Re, 
VS DT OT A Re, (15) 


7 eens 


The friction coefficient and the Stanton number are evaluated by 
‘the basic wall shear and heat flux laws respectively. Both of these laws 
are based on incompressible flow along a smooth, isothermal, impervious, 
flat plate. The friction coefficient and Stanton number are: 


for laminar flow 


Cc 
f,2,0 | 9.245 (16) 
6 
Chot.0 Am 0.22 : , (17) 
pr? Re, 
And for turbulent flow 
Cee — 9.245 , 9:0723 Re, 


_Ls 
“e * Res * TOO F Re, (109 gReg) (18) 
022 22 0.0123 
c eer 0g rey) (19) 
htop» Bre, Pr? (100 + Re 10 


In order to account for the various boundary layer effects, the 
Stanton number and friction coefficient given above are modified by the 
corresponding influence coefficients: 

z 


C = ¢ i | for x = h,f 
Xsy Xs¥s0 Xs¥oZ y= hit (20) 


The influence coefficients are shown by I where the subscripts 


XVoz 
x and y indicate whether the influence coefficient pertains to heat or 

momentum transfer (x = h or f) and laminar or turbulent flow (y = % or t), 
respectively. The subscript z indicates the type of phenomenon for which 


the basic laws are being modified. 


. Four phenomena are considered by MEIT. These phenomena and their 
corresponding z-subscripts are given in Tzeble 1. The influence coefficients 
corresponding to each of these effects are formulated below. 
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TABLE 1 


7 . | Phenomena z-subscript of influence 
_ coefficient 
: acceleration 
transpiration 
boundary layer properties 
roughness 
Acceleration 
In laminar flow, Ip og * (1 : 3p)*3, B > 0 
leg 8 = 1.0, 6<0 
. __ us 
Th 2,8 a + 48)", B > 0 
; —Thyae = 1.0, B<0 
, du - . - 
° = 2€ e 
where: 6 = Ue de 
s 
- 2 
f: J Poliggh ds 


In turbulent flow, only the first terms of the appropriate basic 
turbulent laws are modified by the above influence coefficients. 


' Transpiration a3 
‘ Blowing effects are modeled with adaptations to film theory 
(Reference 7). 


ate 


Te Ta area MAREE Aa ne GOR RE AT AS 


|e eee BEIT PTY SEEPS, EL AF aD ten NR ARBs WE Swe hematite, Ham MENS San ao nea recat ca SEGA rien 3h 
oar Mp Rte Mi? eo mee Ra Sat a ce ower ER ENN Mini A AEC ARCA: Miia bmp iaken 
e IRF Yon A te A be oe OOM Reet eee oe 


gn(1 + 2a, R‘B') 


I ol a SEE i LE y22,t 
h,y,B' Any R'B 2 


onl + 2d_ ) RIRB') y= Bt 
s 
Ay RR | 


Tp y,B" = 


where: R CH = Reynolds Analogy Factor 
Cela. - 
R' = Mass to heat transfer coefficient ratio 
d = Blowing reduction parameter. 
R is a dependent variable which is evaluated during the solution 


process. Both R' and xy are input (see Section 4). MEIT, however, does 
provide build-in default values for A, ¥. These default values are: 


Ane = reg = 0.5 


Boundary Layer Properties 
Boundary layer properties of density, viscosity, and Prandt] number 
are evaluated at the reference enthalpy h’' 


h' «= ah, + bh, + ch, 


The property influence coeffic*ents are: 


8 d e i 
I -.= 2 fut = .- vs tg OS ’ 
I, yep (-) (i) » X Fshs ¥- ay ; . = <f 


where the constants a, b, c, d, and e for all combinations of x and y are 


given in Table 2. 
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‘TABLE 2 
Constants Used to Evaluate Property Influence Coefficients 


Property 
influence 
coefficients 


. 


SHON IRAP oe A eT I 
$ - ‘ . “TERS € 
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Surface Roughness 
In laminal flow, 


: : In.tur : le er ce al 


SOS SERIE The 9 Spe BETAS OT RC 


In turbulent flow, the influence coefficient due to surface rough- 
‘ ness is based on correlations from PANT and Stanford heat transfer data 
(Reference 8 and 9), as well as Stanford and NSWC friction data (Refer- 
ence 9 and 10). The turbulent influence coefficients due to roughness 
are: | 


} 


le ep % 140.5 £(k/e) 9(X) 


| 


lace = 1 + 0.3 (k/@) 9(X) 


where: f(k/e) = 1+ 0.09 (k/e) + 0.53 (1-e7*/®) 
OX) = X+1.5 (eX) for x>0 
= 0 . | for Xx <0 


X = log “T5-5 
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+ = Pellek Ue [Py 


Tee ‘fe {'w t,s 
Ye Ww \Pe 7p 2 


CF ts smooth wall friction coefficient given by equation (19) ~ 


k = surface roughness 


k is input either as function of material or location (see Section 4). 

The above formulations for H, F, Ch» Ce h are only for Taminar and 
turbulent flows. To evaluate these four parameters for transitional flow 
the following relation is used: 


P= (1- f) P, + fP, - (21) 


where P is one of the four parameters above and f is the transitional 
intermittency factor. 


The transitional intermittency employed in MEIT is based on the 
work of Persh (Reference 11), and according to the interpretation of Dahn ii 
(Reference 12). 


a 


ee 
2 
Reg (Cee - Coy) 
2 
where a = REG typ (Cet - Ce eer 


and the subscript tr refers to conditions at the transition point. 
f is set to zero in laminar flow, unity in turbulent flow, and 
varies between 0 and 1 in transitional flow. 


---~- —This completes the formulation of the four required apart used | 
in the solutions of equations (1) and (2). 
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. is generated by the program to ensure adequate solution accuracy of the 
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2.3 SOLUTIONS OF BOUNDARY LAYER INTEGRAL EQUATIONS 


As mentioned above, the required input to MEIT are surface shape, 
boundary layer edge conditions, boundary layer gas properties, and wall 
conditions. These quantities are input in terms of body points. A finer 
grid in terms of integration points, which include all the body points, 


integral equations. The boundary layer edge conditions, gas properties 

and wall conditions at each integration point are obtained by linear inter- 
polation from the input. The solution procedures of the boundary layer 
integral equations (1) and (2) consist of: 


e start-up series solutions at the first three integration points 


e finite difference numerical solutions for the rest of the 
integration points 


The solutions at the first integration point are given by: 


0.245 v7 (1+ R,B ?) z 


e,°3 —HC 
] (3 + H) Ue f haz 
Si 
- ni 22 v, (1 + By ) 
ae 1 


The solutions at the second and third integration point are related 
to the first integration point by: 


@ = (1 + ay?) 


6 =o (1+ by) 


13 +H, 7 (0-659 (3 + Hy) @ (y-1) _ (H)-0.614) a (y-1) 
Ol SS de, 
he ee eee ee ee ik 
where: a sz a 
Ke mh E " 1 2x 0.659 a (y-1) 7 2(1-F, )(1+B') a(y-1) | 
v= S/Rrog 
1 - (p3/p,) 
and oo, 
3 


In tie above formulation, Ruef is an arbitrary constant radius and 


y is the isentropic exponent. The subscripts 1 and 3 denote the first and 
third integration point condition respectively. 


The solutions of the first three integration points were originally 
derived for reentry nosetip ablation prediction, and consider the first 
integration point to be a stagnation point. In MEIT, these solutions serve 
only as the start up procedure and should be ignored in the output. 


For the rest of the integration points, the following implicit 
finite difference scheme is used: 


Pet 7 Fxygay + O-5Fy py t Fe lsy = $y.) eseean 
where: Fe =r Pal, *6 (see Equation (1)) 
FAY Pag (My 4 - ho (see Equation (2)) 
dF 
ee eae 3 
tas 
aoe 
h ds” 
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I is the integration point index and Fe and FA are both evaluated 


from equation (1) and (2) respectively. 
Since the values of Es I depend on FY I the solution is obtained 


by iteration. This iteration is local because closure is obtained at each 
integration point before proceeding down the body to the next integration 
point. Convergence is based on changes of less than 0.1 percent in both 
the heat and momentum transfer coefficients between successive iterations. 
If the iteration fails to converge in 30 tries, a local explicit solution 
is obtained by setting Fit = x,I-1° and subsequently reevaluating Fy, 


based on the resulting value of ro I’ before proceeding to tne next 


integration point. 


SECTION 3 
VERIFICATION OF MEIT 


Three computer progranis were used in the check out of MEIT. ‘These 
three codes are: 


e@ ARGEIBL 
e BLIMP 
e TBL 


Both ARGEIBL and BLIMP have no provisions to account for surface 
roughness. ARGEIBL solves only the et.argy integral equation and the heat 
transfer coefficients it calculates are routinely multiplied by 0.75 to 
improve their accuracy. BLIMP uses an exact implicit technique to solve 
the boundary layer differential conservation equations. The solution from 
BLIMP is believed to be the most accurate of these codes and was used as 
the primary criterion to determine the validitv of MEIT in the check out. 
TBL is developed by Pratt and Whitney Aircraft Company and, like MEIT, it 
uses both the momentum and energy integral equations to calculate the 
Stanton number. It also has build-in routines to account for surface 


roughness. 


The check out is performed by generating and comparing heat transfer 
coefficients for two typical rocket nozzle en.‘ronments using these four 
computer codes. The selected rocket nozzle environments are the Aerojet MX 
upper stage and the C/CAN nozzle. The results of the check out and the 
conclusions and recommandations are summarized below. 


3.1 AEROJET MX UPPER STAGE NOZZLE 


An analysis of the Aerojet MX upper stage nozzle by BLIMP, ARGEIBL, 
and MEIT as used to validate MEIT. The geometry of the nozzle is shown 
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in Figure 1. A listing of the nozzle wall coordinates is given in Table 3. 
The nozzle has a 3.146-inch radius throat aid it uses a PEG/FEFO propellant. 
Tne elemental composition of the propellant and the associated chamber 
conditions are given in Table 4. All this information was obtained from 
Reference 13. The edge conditions were calculated by the ACE code 
(Reference 14), and the results in terms of edge pressure, temperature, 
enthalpy, and velocity are presented in Figure 2. The nozzle radius is 
Ppiotted as a function of stream length in Figure 3. A wall temperature 

of 5500°R was used in a]}1 the following calculations. 


Smooth wall solutions were obtained using MEI7, ARGEIBL, and BLIMP, 
and the heat transfer coefficients from these solutions are compared in 
Figure 4. The MEIT and ARGEIBL solutions are very close to each other for 
regions downstream of the throat. Both of these solutions are high com- 
pared to BLIMP results. As mentioned above, the ARGEIBL solution mus* be 
reduced by a 0.75 factor to accurately predict the heat transfer coeffi- 
cient. Using 0.75 as the modifier, the modified MEIT and ARGEIBL solutions 
are again compared to the BLIMP solution in Figure 5. From this figure, it 
can be seen that both the reduced MEIT and ARGEIBL solutions compare much 
better with the BLIMP solution. This shows that MEIT, like ARGEIBL, re- 
quires a modifier of 0.75. 


For the unblown, smooth wall and turbulent flow calculations (as in 
the above case) it can be shown from Section 1? and Reference 1 that MEIT 
is solving the same Energy Integral Equation as in ARGEIBL (i.e., for this 
case the Momentum Integral Cquation does not enter into the solution pro- 
cess). It is therefore of interest to find out the reasons for the dif- 
ference. in the MEIT and ARGEIBL heat transfer coefficient solutions. The 
chief reason is the different formulation for the local Stanton number in 
the two codes. In MEIT, the local Stanton number is given by equation (19) 


; a 0.22 _ , 9.0123 Re® _ (194 pe,)""-5| 
h,t,0 Pr”? Re, Pr* (100 + Re,) 


and modified by the influence coefficients due to acceleration and boundary 
layer properties. In ARGEIBL, the Stanton number is defined by: 
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TASLE 4. PEG/FEFO PROPELLANT DATA 


ELEMENTAL COMPOSITION 


ELEMENT 
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2.39527 
1.13509 


1279248 


2.27887 
0.07275 
0.68503 
0.07071 


cee a RE AT SS ae ene ae RM en oY 


"alzzou a6e3s uaddn yy yefouay ayr 


30 SUOLJLPUOD aBpzZ °2 aunbLy 
S3HONI HLGNZT PYTSLS 
en . r zs s : z : 
Ftp conrrcinasaline LA = i) 
ans 
, z 
2 
» § 
e 5 
> m 
§ 3 
aoe 
oe 
: : Oo8 
- 
= 
e . e ° 
oie sities ce Lal Wee cea ce ok wrath ne ener Ae alte “Bil mabdeDeadn okD int San an a cbs tent Net Bie Nate oo RE tna ln ne nas te are haat ebb 


O0€- 0 
904 200% 
m 
- 
z : 
00s e  S00'y 
< 
= 
£ 
gs 
os 6 000'9 
5 
& 
o0c's 000'9 
o02"s 000s 
a 


#9398/14 (3M) ALIDONZA 


benek 4 


Sot 


co0'v 


000s 


o00'S 


000s 


We (an) BeRLvY3GH3L 


: 
*a_zzou abeys uaddn yy Jefousay ayz wos Yyrbua, weauzs “SA SNp2 a_zzZoN “Ee aunByy 
S3HONI ‘HLON37 Wv3uLS ; 
et o rh zt OL 8 9 ’ z 0 ; 
0 
z on 


24 


GL¥61-v 
v 


t -) 
S3HONI ‘snowy 


On 


S3HON! ‘HLONS] WVSHLS 


gt bet ot 8 9 y 2 0 


25 


GWAITG amu 0 eee 


O3S2L4/N81 ‘LN31D1SS3OO YS4SNVUL iV3H 


TIShol-v 


*a_zzou abeys uaddn yyw yafouey ayy u0y 3U9}9}33909 uaysuesz Yeay “py aunbyy : : 
a 


Be ph tele WARS se EA aed g ia I N ean ee tine ite AR aah! ine OAR th RN oer RAR ct? Bett ait ike mth Stns Bes cn 


Hhiesi sidaacnnbicinn., We aan eB 


*a_zzou abeqs 
daddn XW Jafouay a4 wos UaPI}44a09 waysuedg WRdaY PILsEPOW “§ aunbhy 


S3HON! "HLON3T WY3ULS 4 
ot v1 zi o1 8 9 7 o,. a ct 
) 4 : 
4 
: ot * 
= i ~ ‘® 
om } oon ‘ 
o faa ‘ 
so ok es 
Oo : ‘ 
x . 2 
$ 1 4 
_ . 
+ 4 
ot > ; 
> 4 “ 
5 mn 
ry 5 
S Qo 
bs re) 
we m 
sh 
Q 
om 
Zz 
xs 
Cc 
o 
dNI1G smeaes 4 
TGIZOUY eee ° . 
13 <—oo ae > 
\- 
\. 
\ a \ x 
s “ XN eae , 


2 UM er TH DET EN AIRES PRET OE ONAN MUR SCO PE ROR AD POT 


PRYOR ORE ERP IIE ne RE EERE HE 


+ Ve 
C * 0.0130 (2 ) 
pnt! + aNeres 
(Pr')% Re™* 
where the subscript ‘ denotes properties evaluated at the reference enthalpy. 


Other reasons are differences in the subroutines built in each code. 


These are the following: 


e@  ARGEIBL uses the trapezoidal rule to calculate the numerical 
integration whereas MEIT uses an averaging technique. 


e ARGEIBL uses a cubic curve fit to interpolate the required 
' properties whereas MEIT uses a linear interpolation. 


To check if MEIT is indeed solving the same energy integral equation 
as in ARGEIBL for the unblown, smooth wall and turbulent flow calculation, 
the Stanton number formulation in MEIT was replaced by the ARGEIBL formula-. 
tion, and a solution was obtained using this new 'ARGEIBL-MEIT' code. The 
heat transfer coefficients from this calculation are again compared to those 
given by ARGEIBL in Figure 6. The two solutions compare very well in the 


" mose cap region and the differences in other locations are basically due to 


the different. inherent routines mentioned above. 


In order to check the surface roughness formulation, a rough wall 
solution was generated by MEIT. A roughness height of 2 mil was used 
throughout the nozzle. The roughness augmentation from this calculation . 
is compared to those given in Reference 15 in Figure 7. The results from 
Reference) J5 are based on local energy thickness and conditions given by the 
smooth wall BLIMP solution and were hand-calculated using the MEIT surface 
roughness formulation. Figure 7 shows that these two solutions compare very 
well with each other. The higher values given by Reference 15 are mainly 
due to the use of smooth wall momentum thicknesses in the hand calculation. 
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| 3.2  C/CAN NOZZLE 


The C/CAN nozzle is used primarily to compare the rough wall solu- 
tions generated by McIT and TBL. The geometry of the nozzle is shown in 
Figure 8. A listing of the wall coordinates of the nozzle is given in 
Table 5. The nozzle has a throat radius of 0.960 inch and uses a HT 90/18 
propellant. The elemental composition of the propellant and the chamber 
—— conditions are given in Table 6. The above information was obtained from 
. Reference 16. The edge conditions were again calculated by ACE and are 
presented in Figure 9. The nozzle radius vs. stream length plot is given in 
Figure 10 and the wall temperature is assumed to be 4939°R. 


| ae , Both smooth and rough wall calculations were performed by MEIT 
i Smooth and rough wal? solutions generated by TBL were furnished to Aero- 
therm by Atlantic Research Corporation (Reference 16). The rough wall 
calculations were based on a roughness height of 2 mil. The solutions are 
: : compared in Figure 11. For the smooth wall case, the MEIT and TBL solutions 
is are aimost identical. For the rough wall case, MEIT predicts lower heat 
7 : transfer coefficients than TBL except for locations just downstream of the 
i combustion chamber. According to Reference 16, the TBL solutions, like 
1° those generated by ARGEIBL, are normally modified by modifying factors to 
make them compatible with measured values. These modifying factors are 
ua developed from actual motor firing data. However, since there were no fir- 
: ing data available at the time of the check out, the question of whether 
MEIT or TBL is more accurate in rough wall prediction cannot be answered 
at this time. 
Reference 16 aiso indicates that the roughness modelling in TBL 
Nee <4 would not approach the smooth wall modelling when the surface roughness is 
a reduced. To make sure that this is not the case for MEIT, an additional 
MEIT calculation was performed using a roughness height of 0.0] mil. The 
results are exactly the same as the smooth wall solutions shown in Figure 
11. This indicates that MEIT is self-consistent. 


3.3 CONCLUSIONS AND RECO}MENDATIONS 


The findings and conclusions based on the above discussion are 
summarized below: 
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TABLE 6. HT 90/18 PROPELLANT DATA 


Pchamber = 68.027 atm 


"chamber = 6627.6 °R 


ELEMENTAL COMPOSITION: 
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For nonblown , smooth wall and turbulent flow, both MEIT and 
ARGEIBL solve the same energy integral equation. 

Using the BLIMP solution as the criterion and a modifying fac- 
tor of 0.75, the modified MEIT solution is better than the 
modified ARGEIBL solution. 


MEIT and TBL predict essentially the same heat transfer coeffi- 
cient for smooth wall calculation. — 

For rough wall calculation, MEIT predicts a lower heat transfer 
coefficient than TSL except in regions inmediately downstream 
of the combustion chamber. : 


MEIT is self-consistent in its roughness modelling. 


The recommendations are: 
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For both smooth and rough wal] calculations, a factor of 0.75 
should be used to modify the heat transfer coefficient generated 
by MEIT. 

Experiments or analyses should be conducted to deduce the heat 
transfer coefficient from actual motor firing data to check the 
validity of the surface roughness modelling in MEIT. 


A technique needs to be devised to evaluate nozzle material 
roughnesses for input into MEIT. 


When there {fs no firing data to modify the TBL heat transfer 
coefficient, the MEIT solution modofied by a 0.75 factor should 
be used for rough wall predictions. 
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SECTION 4 
DESCRIPTIONS OF INPUT AND OUTPUT 


This section provides detailed user oriented input instructions and 
a description of the output. The input instructions are presented in 
Section 4.1 and the output features are covered in Section 4.2. 


4.1 INPUT INSTRUCTIONS 


The program input consists of three card sets: 
1. General program constants and transition information 
2. Surface shape, boundary layer edge and wall conditions 
3. Boundary layer gas properties 


The description of these three card sets is given below. 


Card Set 1— Control Cards 


This card set is made up of three cards which provide general pro- 
gram constants and transition information. The four blowing reduction 
parameters, Ax,y? which are used to calculate the blowing influence coef- 


fients (see Section 2.2), are also included in this card set 


Card 
No. Column Format Data Unit 
J 1 1) —s NSM - No. of material, up to 3 allowed - 
2 TROUGH - Roughness height input flag / - 
0 - Input as function of material 
(see below) 


1 = Input as function of location 
(see Card Set 2 
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i 


| 
| 
| 
| 
: 


[ 
| t 
| E Card 
é. 
| » 
| " 3-10 
—_ 
i 11-20 
a 2 1-10 
t 11-20 
21-30 
a 31-40 
i 
i 41-50 
i 
i 
t 
f 
} 2+ NSM 142 
& e- 
' 2-4 
8-6 
7-8 
* 9-12 
13-22 


subsets are: 


e 2ZsP 


No. Column Format 


F8.5 


F10.5 


F10.5 
F10.5 
F10.5 


F10.5 
F10.5 


Repeat card 2 for each material (number of cards equals NSM). 


12 


12 


12 


12 


@  IMAT - Material {nan default value = 1 
- Axial wall coordinates (in.) 


freak le A Ae Cite ad aah LM a, BAST Rn oe ne TN a en eta Debt EA 


BTS - , > default value = 0.5 ais 
BTH - Ah, t? default value = 0.5 - 
BLS - van default value = 0.35 - 
—~BLH = Ay ye default value = 0.35 - 
RUFL - Surface roughness height as func- mil 


Data | Unit 

CMH - Hass to heat transfer coefficient 3 
ratio, default value = 1.0 

GAM2 - Isentropic exponent, default value 


tion of material, input only if 
IROQUGH = 0 


NS - No. of body points (maximum 60) 


IBRUPT - Abrupt transition flag ~ 
0 - laminar and transitional flow 
1 = fully turbulent: flow. 


NREYCR ~ Transition flag - - 
»7 - all turbulent: 
1 = all laminar 
4 = transitional 
IPRNT - Output print flag. = 
1 - Detailed output at body point 
2 - Detailed output at integration 
point i: 
Blank | 


DLTRAN ~ Axial location at which transition 
takes place, input only if NREYCR = 4 in 


Card Set 2 — Body Point Data 


This card set is made up of nine subsets and they provide the sur- 
face sharp, edge and wall condition information to the program. The nine 


Soetetaig Kintiy 22 pt hen he ges i ican aca Na EO Bien a le nae eect me we we ‘evstiabe catiaocll 


petagal werd tpetoe 


e RSP = radial wall coordinates (in) 

e PE ~- Edge pressure (atm) 

e@ HE - Edge enthalpy (B8tu/1bm) 

e UE - Edge velocity (ft/sec) 

e TSP - Wall temperature (°R) 

e BPSP = Blowing parameter = Wo u Gy 


@ RUF2 - Surface roughness height (mf1) 
'  RUF2 is input only if IROUGH = 1 (see Card Set 2). 


The nine subsets are input in the above order. Each of these sub- 
sets has NS entries and each entry corresponds to a point on the heated 
surface of the body. The input format of all the subsets are F10.3 except 
IMAT which uses an input format of I10. The values of each subset are 
entered eight to a card. Since the solution for the first body point is 
based on stagnation conditions (see Section 2), it is recommended that the 
first entry in the velocity subset (UE) be zero. 


Card Set_3 — Gas Property Tables . 

The gas property table provides the thermodynamic and transport 
properties of the boundary layer gas. This table is divided into sub- 
tables based on pressure. Within each pressure subtable, they are ordered 
either on temperature or enthalpy. If they are ordered on temperature, 
the same temperature array must be used for all other pressure subtables. 
The same is true for the enthalpy entries {if the properties are input as 
functions of enthalpy and pressure. 


Card 
. No. Column Format Data Unit 
4 7 — 2 T. Te3 13 ___IPMAX. = No. of pressure entries (maximum 21) 
| 4-6 13 IHMAX - No. of temperature or enthalpy 
entries (maximum 35) 
7-12 13 LTB - Temperature enthalpy flag 
Q- Properties input.as functions of 
pressure and enthalpy 
T- Properties input as functions of 
pressure and temperature 
40 
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No Column Format Data Unit 
2 1-10 10.3 PT -< Pressure atm 
_ 3 to TFMAX + 2 _ - | 
1-12 E12.4 HT -  Enthalpy -Btu/1bm 

13-24 £12.4 TEMT - Temperature oR 
25-36 £12.4 EMT - Molecular weight 
37-60 Blank ; 
61-72 E12.4 ET - Viscosity lbm/ft-sec . 


. 73-80 F8.4 GT - Prandtl No. 


Cards 3 to IHMAX + 2 are repeated until the total number of pressures 
is equal to IPMAX. Note the pressure, enthalpy, and temperature entries must 
be input in ascending order. 


4.2 OUTPUT DESCRIPTION 


MEIT output can be divided into three categories. These are: 
e@ Output of Input | | 
® Output of Calculation Results 

e@ Debug Output 


The descriptions of each of these output are given below. 


Qutput of Input 


The program output begins with the output of the input. This output 
is made up of three parts: i . 


e General program information 
e Thermodynamic table 
e General input information table 


The general program information prints out the program constants 
contained in the first three input cards. These include mass to heat trans- 
fer coefficient ratio, isentropic exponent, the four blowing reduction 
parameters, and transition information. The FORTRAN names corresponding to 
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these variables are given in Section 4.1. The thermodynamic table tabulates 
the thermoaynamic and transient properties of the boundary layer gas as 
functions of oressure. This table is self-explanatory. The general input 
information table gives the surface shape, boundary layer edge conditions 

and wall conditions. The variables printed in this table, their FORTRAN 
names and definitions are listed below: 


— BODY PT NO (J) 
Index of the body points 

-: INiEG PT NO (1) 
Index of the integration point for which the computed parameters 
are printed 


~ MATL NO (MATL) 
Material index. 


— STREAM LENGTH, inch (S) 
Stream length along the heated surface from the initial point 
to the integration points 


— AXIAL LENGTH, inch (Z) 
Axial coordinate of the integration points with respect to the 
initial point 
— | RADIAL LENGTH, inch (R) 
Radial coordinate of the integration points with respect to the 
centerline 
— BODY ANGLE, deg, ees (THETB) 
Angle which the tangent to the surface makes with respect to the 
centerline 
~ NORMALIZED ABLATION RATE (BPSP) 
Normalized ablation rate B' at the wall 
gt = (ev), 
Palaly 


—- WALL TEMP, °R (TW) 
Temperature of the wall 
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PRESSURE, atm (PE) 


Edge pressure of the boundary layer 
ENTHALPY, Stu/Tbm (7) 
Edge enthalpy of th: boundary layer 


VeLOCITY, ft/sec (‘'%) 
Edge velocity of the boundary layer 


SURFACE ROUGHNESS, mil (RUFMIL) 
Surface roughness height of the wall material 


Output of Calculation Results 
The results of the code calculations are printed out in three tables. 
These three tables are: 


1. 
2. 
3. 


Viscous Flow - Edge Properties 
Viscous Flow - Wall and B. L. Recovery Properties 


Viscous Flow - Boundary Layer Solution 


The variables printed in each table, their FORTRAN names and defini- 
tions are listed below. 


1. 


Viscous Flow - Edge Properties Table 


BODY PT NO (J) 


Index of the body points 


INTEG PT NO (I) 


are printed 


STREAM LENGTH, inch (S) 
_Stream length along the heated surface from the initial Pole: to 
the integration points 


VELOCITY, ft/sec (UE) 


Velocity at the edge of the boundary layer | 


MACH NO (HCAM) 


Mach number at the edge of the poate layer 
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_ Index of the integration point for which the computed parameters 


2. 


cata se aaa 


ENTHALPY, Btu/1bm (HE) eo 
Ehthalpy at the edge of the boundary layer - 


TEMPERATURE, °R (TE) 
Temperature at the edge of the kamen — 


DENSITY, Ibm/ft® (ROE) 
Density at the edge of the boundary layer 


VISCOSITY, lbm/ft-sec (VISE) 
Viscosity at the edge of the boundary layer - 


UNIT RE NO, 1/ft (URE) 
Unit Reynolds number at the edge of the boundary layer 


Note that the Mach number given in this table is valid only if 
there are no condensed species in the freestream. 


Viscous Flow - Wall and B.L. Recovery Properties Table 

BODY PT NO (J) 

Index of the body points 

INTEG PT NO (L) 

Index of the integration point for which the computer parameters 
are printed 

STREAM LENGTH, inch (S) 

Stream length along the heated surface from the initial point to 
the integration points 


WALL TEMPERATURE, °R (TW) 
Temperature at the wall 


. WALL ENTHALPY, Btu/lbm (HW) 


Enthalpy of gas at wall temperature and pressure 


WALL DENSITY, Ibm/ft? (ROW) 
Density of gas at wal] temperature and pressure 


WALL VISCOSITY, lbm/ft-sec (VISW) 
Viscosity of gas at wall temperature and pressure 
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‘Momentum thickness 6 of the boundary layer 


_ ENERGY THICKNESS, mil (PHI) 
’ Energy thickness $ of the boundary layer 


_ MOM THICK RE NO (RETH) 
. Reynolds number based on the momentum thickness, Re, = 


RECOVERY ENTHALPY, Btu/1tbm (HR) 


_ Recovery enthalpy defined by h. = he + F (he - he) 


RECOVERY FACTOR (RECOV) V2 
gone. for laminar flow 
Recovery Factor defined by F Pr ,n ={173 ar:.turbuleat lon 


SENSBL CONV HEAT FLUX, Btu/ft? sec 
Sensible convective heat flux defined by ay = PelleCy (h - hy) 


CF/2 
Function coefficient C./2 = t/(p,u,7) 


Viscous Flow — Boundary Layer Solution Table 
BODY PT NO (J) 
Index of the body points 


INTEG PT NO (I) 
Index of the integration point for which the computed parameters 


are printed 


_ STREAM LENGTH, inch (S) 


Stream length along the heated surface from the stagnation point 
to the integration points 


MOMENTUM THICKNESS, mil (THE) 


SHAPE FACTOR (HSF) 
Boundary layer shape factor, H = 6*/6 


where 6* = f .) ( “eee) dy is the boundary layer displacement thickness. 
=: ee eeGre 
6 
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— ENERGY THICK RE NO (REPH) 0.u.6 
Reynolds number based on the energy thickness, Re, = ; £ 
e 


oe HEAT TRANS COEFFICIENT, tbm/ft?sec (RUCH) 
Heat transfer coefficient, Pag Ch, = = 4 /(h, -h ) 


where C, is the Stanton number and aw is the wall heat flux 


~ REYNOLDS ANAL FAC (RAF) 
Reynolds analogy factor C,/(C_/2) 


— _ INTERMITTENCY (ADML) 
Boundary layer intermittency factor f, where 0 < f < 1 for flow 


ranging from fully laminar to fully turbulent regime 


Heat transfer augmentation due to surface roughness, this 
quantity is equivalent to the roughness influence coefficient 


As mentioned in Section 2, the solutions for the first three integration 
points are only for start up purpose. These solutions should be {gnored - 


in the output. 


Debug Output 
There are four error messages output by MEIT. The error messages 
and their meanings are given below. 
= talalalala Wrong Pressure Input Ordering 
The pressure entries in the boundary layer gas properties table 
are not in ascending order. 


t 
i 
j 
{ 
j 
. 
f : _ HEAT TRANS AUGMENT (RUFSMT) 
f 


i 
i 
! 


—  tititieee Wrong Temperature/Enthalpy Input Ordering 
The temperature or enthalpy entries in the boundary layer gas 
properties table are not in ascending order. 


— ‘*eeweeee® Input Temperature/Enthalpy Error 
The corresponding temperature or enthalpy entries in the boundary 
layer gas properties table are not the same for different pres- 
sure entries. - 
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— Computation of MEIT equations did not converge at Point I. 
| Pua The MEIT calculations at Integration Point I did not converge 
| after 30 tries. 


There are also other debug output built in MEIT. However, since 
these output are only useful to the very sophisticated users who are inti- 
mately familiar with the program logic, they are not described here. 
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SECTION 5 - | 
SAMPLE PROBLEMS. 


Presented in this section are two sample problems which were run 
on a Univac 1108 Digital Computer. For each sample preblem, the following 
is presented: 


e A brief description of the nature of the prvblem and solution. 
@ A listing of the input data deck. : 
e A listing of the output. 


Sample Problem 1 e 
This is a smooth wall calculation for the Aerojet MX upper stage i 
nozzle with a PEG/FEFO propellant. The details of this system were 
described in Section 3. A total of 21 body points were input. The wall na 
temperature was assumed to be 5500°R. A IPRNT = 3 option was used which 
printed out only the body point solutions. . \ 
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Sample Problem 2 


A rough wall calculation for the C/CAN nozzle with a HT 90/18 propellant 
was performed. The details of this system were given in Section 3. Fifty- 
six body points were input. The surface roughness and the wall temperature 
used in this calculation are 2 mil and 4939°R respectively. Solutions at 
all the integration points were printed (1PRNT=2). 
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Section 6. Program Listing 


A listing of the Fortran source deck is given in this section. The common 
block and block data are presented first, the driver routine which is the 
main program is given next, followed by all the subroutines in the program. 
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